Dopant-induced crossover from ID to 3D charge transport in conjugated polymers 
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The interplay between inter- and intra-chain charge transport in bulk polythiophene in the hop- 
ping regime has been clarified by studying the conductivity a as a function of frequency cu/2n (up 
to 3 THz), temperature T and doping level c. We present a model which quantitatively explains the 
observed crossover from quasi-one-dimensional transport to three-dimensional hopping conduction 
with increasing doping level. At high frequencies the conductivity is dominated by charge transport 
on one-dimensional conducting chains. 
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The charge transport mechanisms in conjugated poly- 
mers, although extensively studied over the last two 
decades, are still far from completely understood. Not 
only the behavior around the insulator-to-metal transi- 
tion (IMT), which can be induced in several polymer ma- 
terials upon appropriate doping, but also the nature of 
hopping transport in the deeply insulating regime are not 
yet resolved. While some studies indicate that transport 
is dominated by hops between three-dimensional (3D), 
well conducting regions |l],|| , in other cases the strongly 
one-dimensional (ID) character of the polymer systems 
appears to be a crucial factor . 

In investigating the nature of hopping transport in con- 
jugated polymers, studying the temperature and doping 
level dependence of the DC conductivity is an impor- 
tant tool. Since the DC conductivity is determined by 
the weakest links in the conducting path spanning the 
sample, the study of a DC (T) gives insight in the slowest 
relevant transport processes in the system. 

On the insulating side of the IMT, the DC conductiv- 
ity is predicted by many models to follow the well-known 
hopping expression 
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where the value of 7 and the interpretation of T depend 
on the details of the model. The original Mott theory 
for 3D variable range hopping with a constant density of 
states (DOS) at the Fermi energy predicts 7 = 1/4 
while several modifications of the model have been pro- 
posed to describe the frequently observed value 7 = 1/2. 
Studying the dependence of 7 and T on doping level 
c provides the opportunity to discriminate between the 
various hopping models and extract parameters deter- 
mining the conductive properties like the DOS and the 
localization length. 

While the DC conductivity is sensitive to the slowest 
transport processes, the AC conductivity o~(uj) provides 



information about processes occurring at time scales 
t k lo^ 1 . Especially in conjugated polymers, where 
intra-chain and inter-chain transition rates may differ by 
orders of magnitude, knowledge of (j{oj) at high frequen- 
cies can help to clarify the properties of charge transport 
on a polymer chain. 

In this Letter, we present a systematic study of the 
charge transport in a conjugated polymer far away from 
the IMT, as a function of frequency, temperature and 
doping level. By selecting a polymer system with very 
low inter-chain mobility, a separation of inter-chain and 
intra-chain contributions to the conductivity can be made 
when the applied frequency is varied over 12 decades. 
At low frequencies, transport between chains is studied, 
while at frequencies well above the inter-chain transition 
rate, intra-chain conduction is probed. 

Experimental - The experiments were performed on 
discs (thickness between 0.4 and 1.0 mm) of pressed 
powders of the conjugated polymer poly(3,4-di-[(R,S)-2- 
methylbutoxy] thiophene), abbreviated as PMBTh, the 
synthesis of which is described elsewhere [Q . The sam- 
ples were doped with FeCl3 in a dichloromethane solution 
at doping levels 0.01 < c < 0.22; here c is the number 
of doped carriers per thiophene ring, which was deter- 
mined with the aid of Mdssbauer spectroscopy . Af- 
ter doping, the solvent was evaporated and the resulting 
powders were vacuum dried overnight. The conductivity 
of the samples remained unchanged in an ambient atmo- 
sphere for several weeks. Contact resistances were less 
than 5% of the sample resistance |l(]]. The conductiv- 
ity data were taken in the range 5 Hz-3 THz with the 
aid of several experimental methods, which are described 
elsewhere p| . The conductivity in the range 0.3-3 THz 
was determined from the transmission measured with a 
Bruker FTIR spectrometer. 

DC conductivity - The temperature-dependent static 
conductivity of samples with doping levels ranging from 
0.01-0.22 is plotted in Fig. |]. Here log a DC is plotted vs. 
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T -1 / 2 , so that data sets following Eq. @) with 7 = 1/2 
fall on a straight line. In the inset, the value of 7 is 
plotted vs. c. The exponent 7 has been determined by 
plotting the logarithm of the reduced activation energy 
W = d In ajd In T versus InT and fitting it to a straight 
line; the slope of this line gives 7 p2"l| . The data show 
a clear transition in the DC behavior as a function of 
doping level around c = 0.12. For low doping c < c , 7 
values are grouped around 1/2, whereas for c > c , 7 is 
close to the Mott value 1/4. The conductivity data for 
c < c are now fitted to Eq. (jl]) with 7 fixed at 1/2 (solid 
lines), while the data for c > c are fitted with a fixed 
7 = 1/4 (dashed lines). Note that as the conductivity is 
many of orders of magnitude below the minimum metal- 
lic conductivity ~ 100 S/cm and depends very strongly 
on T, PMBTh is far away from the IMT at all doping 
levels. 

Many authors have reported the conductivity in 
conjugated polymers to follow Eq. (|l|) with 7 = 1/2, and 
several models have been proposed to explain this value. 
In disordered systems, the single particle DOS around 
the Fermi energy has a parabolic shape when long-range 
Coulomb interactions between charge carriers are dom- 
inant ]l3| ]. Inserting a quadratic DOS in the original 
Mott argument yields the exponent 7 = 1/2. Data on 
various conjugated polymers close to the IMT have been 
interpreted within this Efros and Shklovskii Coulomb- 
gap model JlJ]. Alternatively, it has been argued that 
polymer materials can be viewed as 3D granular metal- 
lic systems when the strong, inhomogeneous disorder in 
polymer materials leads to the formation of well conduct- 
ing 3D regions separated by poorly conducting barriers 
jl],^| . For granular metals, a DC (T), though still not com- 
pletely understood, is widely accepted to follow Eq. (|l|) 
with 7 = 1/2 @. Close to the IMT both models pre- 
dict a crossover from 7 = 1/2 to 7 = 1/4. Such transi- 
tions have indeed been observed in a number of systems 
0^0. However, because PMBTh is far away from 
the IMT at all c, these models cannot be applied here. 

To explain our data we will use another approach, 
which is an extension of the quasi- ID hopping model 
of Nakhemdov et al. |l8|. In their model, the charge 
carriers are supposed to be strongly localized on single 
chains (or ID bundles of chains). Variable-range hop- 
ping is possible along the chains, while perpendicular 
to the chains with vanishingly small inter-chain overlap 
only nearest-neighbor hopping is allowed. Although the 
quasi-lD model predicts 7 = 1/2 strictly speaking only 
for the anisotropic effective conductivity perpendicular 
to the chains, it has been successfully applied to ran- 
domly oriented polymer systems ||. In a closely related 
approach, the polymer system is viewed as a fractal struc- 
ture with a fractal dimension df slightly greater than one 
||l9t . The static conductivity on such a nearly ID fractal 
is also calculated to follow Eq. ([!]) with 7 = 1/2. 

Model - To include a transition from quasi-lD hopping 
to 3D hopping as a function of doping level, the models 
mentioned above need to be extended. In the quasi-lD 



model such a transition is expected when the transverse 
overlap is sufficiently increased. In the nearly ID fractal 
model, an increase in the fractal dimension would eventu- 
ally lead to a decrease of 7 down to 1/4. So in both cases, 
a transition to 3D hopping is induced by an increase of 
the inter-chain connectivity. Before we present our calcu- 
lations, let us show how this can be qualitatively under- 
stood. Upon chemically doping a conjugated polymer, 
not only charge carriers, but also dopant counter-ions 
are added to the system. The counter-ions locally de- 
crease the inter-chain potential seen by the carriers, and 
thus considerably enhances the hopping rate ||. Since 
the Fermi energy is shifted upwards upon doping, the 
counter-ions must even be considered as hopping sites 
when the distance of the dopant site energy to the Fermi 
level becomes of the order of the hopping energy. In this 
view, the sharply defined crossover at c can be inter- 
preted as the point where the dopant sites start to play 
an active role in the hopping process, thereby enhancing 
the density of hopping sites and making variable range 
hopping in the direction perpendicular to the chains pos- 
sible. 

We now show quantitatively that the conductivity 
in both doping regimes and the crossover with doping 
level can be explained in a single variable-range hop- 
ping model. We assume hops within an energy in- 
terval E over a distance (X, Z) in the combined par- 
allel (a;) and orthogonal (z) directions, with localiza- 
tion lengths L x and L z respectively. For an electron 
on a chain, the local DOS n contains two contribu- 
tions, n from the chain itself and n-^ from the neigh- 
boring chains and from intermediate dopant sites. As 
the latter contribution depends stronger on c than the 
former, ni/rio rises with doping level. Following the 
usual variable range hopping arguments, the conductiv- 
ity oc exp[— 2X/L X — 2Z/L Z — E/k B T] should be maxi- 
mized under the condition 2X(2Z) 2 (n)E w 1, where (n) 
is the DOS averaged over the volume (2X,2Z,2Z), for 
which we write (n) = n (L z /2Z) 2 + n x p3]. We intro- 
duce £ = 2X/L X , C = 2Z/L Z and e = E/k B T and note 
that (n) does not depend on X, leading to e = £, i.e. 
<j oc exp[— 2£ — £]. Optimizing in the 'high doping' limit 
C 2 > njn-y we find 2£ = C and = 4k B TL x L 2 z n 1: 
which gives a oc exp[-4£] = exp[-(T ' ll9 ' l /T) 1 / 4 ] with 
rphigh _ (34/ (^£^£2^^ p or ^ e l i ow _cioping' regime 
C 2 < %/%; we get £~ 2 = k B TL x L 2 z n Q and ( = 0, 
which leads to er oc exp[-2£] = exp[-(T ( 5°" , /T) 1 / 2 ] with 
T iow = 4/(k B L x L 2 z n ). Note that although in this limit 
the T-dependence (7 = 1/2) is determined by the dom- 
inating hops in the chain direction, occasional hops be- 
tween chains will still happen. In the high doping regime, 
hops in all directions are equally likely and become long- 
ranged at low T . 

AC conductivity - In Fig. a(u>) at room temper- 
ature is plotted for three samples with doping levels 
0.03 < c < 0.22. At low frequencies, the conductiv- 
ity is seen to be independent of frequency and equal to 
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the DC value. At the onset frequency ui Q /2ir w 1 MHz, 
the conductivity starts to rise, following an approximate 
power law a ~ lo s with s < 1. An extra upturn in 
the conductivity is observed around a second frequency 
oj 1 /2tt w 10 GHz. The temperature dependence of 
the high frequency (200-600 GHz) conductivity was also 
measured between 4 and 300 K (not shown), revealing 
that the frequency dependence a oc u> s with s w 1.6 
is independent of temperature. The absolute value of 
the conductivity only shows a weak (30%) decrease go- 
ing from 300 to 150 K, and is constant with temperature 
below 150 K. 

As was discussed above, the conductivity of a system 
of coupled polymer chains generally consists of contri- 
butions due to both inter- and intra-chain transport. 
When the chains are only weakly coupled, i.e. the inter- 
chain hopping rate Ti nte r is low, the conductivity cr(u>) 
at frequencies u> 3> Winter is dominated by charge trans- 
port processes within single polymer chains. It is widely 
known that in a ID chain, any impurity causes the states 
to be localized, so the chain has zero conductivity in the 
limit uj -> 0, T -> At T -> and uj > 0, the con- 
ductivity of the chain is finite, stemming from resonant 
photon induced transitions between the localized states, 
and is at low frequencies wr < 1 given by Q 

= ^ e2y F T ^ 2 ln2 (!M) (2) 

where v F is the Fermi velocity on the chain, r is the 
backward scattering time and b is the inter-chain sep- 
aration. Following Ref. p3j, this may be written as 
ct (oj) = (it /2b 2 )e 2 g^Llhuj^ln 2 (1 / 'lot) , where g = n Q L 2 z 
is the on-chain DOS per unit length. At finite temper- 
atures, an extra contribution is present due to phonon 
assisted hopping within the chain. The phonon assisted 
conductivity is given by the ID pair approximation S 

7T 3 

<r x (a;,T) = j^e 2 g 2 Llk B Toj\ii 2 (is ph /Lu) (3) 

valid for frequencies uj below the phonon 'attempt' fre- 
quency v p h- The total conductivity at temperature T is 
now given by the sum of the two contributions, a(ui,T) = 
a a (uj) + cr x {w,T) given by Eqs. (||) and (0). The data at 
ui/2tt > 10 MHz have been fitted with this <r(u>), as is 
shown in the inset of Fig. |^. Here, the dielectric loss func- 
tion e"(u>) — (j(u))/e u) is plotted at frequencies between 
10 MHz and 3 THz, together with the fitting line. The 
fits are excellent at high frequencies lu/2tt > 200 MHz. 
The deviations below 200 MHz indicate that either mul- 
tiple intra-chain hopping or inter-chain transitions have 
significant contributions in the MHz regime, which is con- 
sistent with the variable range hopping description at low 
frequencies. 

Parameter values - From the fits of the AC data, the 
parameters r = 10~ 14 s and v p h = 2 x 10 12 s _1 are 
extracted. While the phonon frequency v p h is in good 
agreement with commonly suggested values of 10 12 s _1 



, the scattering time r is an order of magnitude longer 
than reported values for other (highly conducting) con- 
jugated polymers [ p5p^ ], which is likely due to a smaller 
v F resulting from the low band filling. The typical time 
scales now follow directly from the fits, since phonon me- 
diated hops between two sites within a chain occur at 
rates up to T p h — v p h exp[— 2X/L X ]. With X/L x > 1, we 
have Tph^max ~ 10 11 s _1 , equivalent to a local diffusion 
constant D = 10~ 7 m 2 /s. 

Assuming the localization lengths L x along the chain 
and L z perpendicular to the chain to be independent 
of doping level, they can be determined from the com- 
bined DC and high frequency conductivity data. From 
the samples in the low-doping regime, L x and g can 
be extracted using T^ ow = A/(k B L x g ) w 10 5 K and 
Eq. (||). This gives L x = 10 A, indicating that carri- 
ers are localized on the chains in regions consisting of 
two to three rings; furthermore g — 0.1 levels/(eV ring) 
for c = 0.03. For the typical hopping distance along 
the chain we find X = {L x /A){T^ ow /T) 1 ' 2 = 50 A. At 
the onset of the high-doping regime C, 2 ~ 2n /n 1 ; using 
Z w b and T^ gh = 64/(fc B L a; L 2 n 1 ) » 10 8 K, wc find 
L z = 1.4 A, close to reported values for other conjugated 
polymers lying deeply in the insulating regime ]5|,^7j , and 
n 1 = 2 x 10 26 cV _1 m -3 , implying a DOS of 0.6 states per 
eV per dopant molecule. Within the localized regions, 
the carriers move in the chain direction with velocity v x , 
which can be estimated using the fact that in ID con- 
ductors the mean free path l x — v x t ~ L x ; this gives 
v x ~ 10 5 m/s, similar to values observed in ID-organic 
conductors |2^] and other conjugated polymers p9[ . 

In summary, we measured <t(uj) in a conjugated poly- 
thiophene with small inter-chain overlap. We developed 
a model that allows a consistent analysis of the criui) data 
in terms of inter- and intra-chain transport. From the low 
frequency results we have found that carriers are strongly 
localized on ID chains with L z = 1.0 A, and no 3D metal- 
lic islands are present. The high frequency data show ID 
transport along the polymer chains with a scattering time 
t = 10~ 14 s, while intra-chain phonon assisted hopping 
proceeds at rates T p h < 10 11 s _1 . 
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FIG. 1. DC conductivity as a function of tempera- 
ture for several doping levels. The data are fitted to 
°dc = a o ex P[~ (lo/T) 7 ]. For the drawn lines 7 = 1/2, 
whereas for the dashed lines 7 = 1/4. The inset shows the 
c-dependence of 7 and defines the symbols used in this graph. 
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FIG. 2. g'(uj) between 5 Hz and 3 THz at 300 K for three 
doping levels 0.03 < c < 0.22. In the inset the loss func- 
tion e"(uj) — a'(io)/e io is plotted for the same three samples. 
Drawn lines are fits to Eqs. (^|) and (|^). Data are corrected 
for the relatively small dipolar response of the alkoxy side 
chain (the maximal value of e" = 0.1 is centered around 0.5 
GHz), most clearly visible in the undoped polymer. 
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